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ABSTRACT 


LASA  data  has  been  analyzed  to  determine  the  probabilities  that  M  -  m.  and 

o  b 

Modified  Spectral  Ratio  (MSR)  can  be  effectively  applied  to  seismic  events.  The  proba¬ 
bilities  have  been  estimated  as  a  function  of  magnitude  and  whether  the  event  is  an 
explosion  or  an  earthquake.  In  the  case  of  the  earthquakes,  the  probabilities  go  from 
about  1. 0  to  0. 0  as  goes  from  5.  0  to  4.  0.  The  drop  is  slightly  faster  for  Mg  -  m^. 
MSR  behaves  similarly  for  explosions.  However,  the  probability  for  Mg  -  in  the 
case  of  explosions  drops  to  zero  roughly  over  the  range  of  from  5.2  to  4.5.  In 
addition,  a  limited  study  of  joint  properties  of  Mg  -  and  MSR  has  indicated  that  they 
operate  quite  independently  in  the  sense  that  if  one  of  the  discriminants  yields  no  decision 
concerning  a  particular  event,  the  probability  that  the  other  can  make  a  decision  is  not 
significantly  affected. 
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LASA  DECISION  PROBABILITIES  FOR  Mn  -  itl  AND  MODIFIED  SPECTRAL  RATIO 

S  D 

I.  Introduction 

It  is  now  generally  recognized  that  if  an  underground  nuclear  explosion  and  a 
normal  depth  earthquake  radiate  equal  amounts  of  energy  as  short  period  compress  - 
ional  body  phases  then  the  earthquake  will  normally  generate  significantly  more  long 
period  surface  wave  energy  than  the  explosion.  One  crude  measure  of  the  surface 
wave  energy  is  the  surface  wave  magnitude,  M  .  The  body  wave  magnitude,  m  , 

O  JJ 

is  a  measure  of  energy  in  the  compressional  body  phases  (1),  (2).  If  M  is  plotted 

against  for  a  set  of  explosions  and  earthquakes,  the  points  tend  to  separate  into 

two  groups.  A  new  event  can  be  identified  as  an  explosion  or  earthquake  by  noting 

where  it  is  located  in  the  ,  nr  plane.  This  is  the  M0  -m,  method  for  discrimi- 

S  b  S  b 

nation.  A  number  of  researchers  have  reported  upon  the  efficacy  of  this  discriminant 
(3),  (4),  (5),  (6),  (7),  (8). 

Most  studies  of  M„  -  m,  have  emphasized  its  value  as  a  discriminant  when  both 
S  b 

Mg  and  can  be  measured.  It  is  also  important  to  determine  the  probabilities 
that  the  appropriate  measurements  can  in  fact  be  made.  In  particular,  if  an 
measurement  can  be  made  at  LASA,  what  is  the  probability  that  the  Mg  -m^  criter¬ 
ion  can  be  applied  using  LASA  data  only?  The  probability  will  depend  upon  the 

value  and  true  nature  of  the  event.  Some  estimates  of  such  probabilities  have 
been  obtained  for  LASA  and  are  presented  in  the  sequel. 

Another  discriminant  of  particular  value  which  uses  LASA  data  is  spectral 
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ratio  (9),  (10).  One  basic  form  of  this  discriminant  makes  use  of  the  ratio 

1. 95 

/  I  X(f)|  df 

1. 45 

P  =  _ 

0.  85 

/  I  X  (f)  |  df 

0.  35 

where  X  (f)  is  the  Fourier  transform  of  10  seconds  of  short  period  data  which  starts 
just  before  the  onset  of  the  P  phase.  If  p  is  plotted  against  m^  for  explosions  and 
earthquakes  the  points  for  the  two  populations  tend  to  separate.  However,  it  is  clear 
that  p  can  be  significantly  affected  by  background  noise.  For  this  reason,  Lacoss  (8) 
introduced  a  modification  to  this  basic  spectral  ratio  test.  The  modification  essen¬ 
tially  consists  of  making  no  decision  on  events  which  have  insufficient  signal  to  noise 
ratio  in  the  two  frequency  bands  used  to  compute  p  .  The  modification  effectively 
eliminates  the  possibility  of  erroneous  classification  of  events  on  the  basis  of  noise 
rather  than  signal. 

The  modified  spectral  ratio  shares  an  important  property  with  the  Mg  -  m^ 
criterion.  If  no  surface  waves  are  visible  on  long  period  seismograms,  one  does 
not  blindly  measure  noise  and  thus  attribute  a  false  Mg  value  to  an  event.  In  a  simi¬ 
lar  way,  the  modified  spectral  ratio  (MSR)  prevents  one  from  blindly  measuring 
p  for  noise  dominated  events.  Thus,  as  is  the  case  for  Mg  -  ,  there  is  a 
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probability  that  MSR  can  be  applied.  Such  probabilities  have  now  been  estimated  for 
MSR. 

Given  any  set  of  events  to  be  classified,  there  is  one  subset,  Cq  ,  to  which 
Mg  -  cannot  be  applied  and  another  subset,  ,  for  which  MSR  indicates  that  no 
decision  should  be  made.  It  is  clearly  of  interest  to  know  the  extent  of  the  intersection 
of  these  two  sets.  Some  data  of  this  type  have  been  gathered  and  are  discussed  in 
this  report. 

The  data  base  for  this  report  is  the  same  as  that  used  previously  by  Lacoss  (8). 
The  explosion  population  consists  of  35  Soviet,  French,  and  American  presumed  or 
known  underground  tests  recorded  teleseismically  at  LASA.  The  earthquake  popu  - 
lation  consists  primarily  of  about  150  events  with  epicenters  in  the  Sino -Soviet  Bloc 
area.  The  material  in  this  report  is  an  extension  of  results  previously  presented  by 
Lacoss  (8). 
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II.  Probability  that  Mg  -  and  MSR  Can  Be  Applied  to  Earthquakes 
In  this  and  all  other  sections  only  events  with  detected  short  period  body  phases 
have  been  considered.  Figure  1  shows  Mg-  m^  data  for  such  a  set  of  events.  A 
possible  decision  line,  which  we  consider  to  extend  linearly  outside  the  range  of  the 
figure,  has  been  drawn  on  the  figure.  The  implied  decision  rule  is  to  decide  explosion 
if  an  event  plots  below  the  decision  line.  Note  that  tills  decision  does  not  require  that 
a  surface  wave  be  recorded  for  the  event.  If  M  can  be  measured  and  the  event  is 
above  the  decision  line,  then  the  earthquake  decision  is  made.  However,  if  noise  or 
other  events  obscure  an  event  and  it  cannot  be  constrained  to  be  in  the  explosion 
region,  no  decision  can  be  made.  Using  the  above  decision  rule,  only  one  earthquake 
in  our  population  would  be  erroneously  identified  as  an  explosion.  The  probability 
that  no  surface  waves  will  be  recorded  for  an  earthquake  and  that  the  noise  or  an 
interfering  event  can  bound  the  event  in  the  explosion  region  appears  to  be  negligible. 
Thus  the  probability  that  M  -  m  can  be  applied  to  an  earthquake  is  the  same  as  the 

u  D 

probability  that  surface  waves  can  be  detected. 

Figure  2  shows  the  basic  data  used  to  estimate  the  probability,  P^m^),  that 
Mg-m^  can  be  applied  to  an  earthquake.  PQ(m^)  is  the  product  of  P^,^(m^),  the 
probability  that  M  can  be  measured  when  there  is  microseismic  noise  but  no  inter- 

J 

fering  event,  and  PjQ(m^),  the  probability  that  no  interfering  event  will  obscure  the 
measurement  of  Mg.  It  is  these  component  probabilities,  P^,Q(m^)  and  P^(m^), 
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which  can  be  directly  estimated  using  the  Figure  2  incremental  histograms.  P^TQ(m^) 
is  estimated  using  the  top  two  histograms.  P^q(M)  is  estimated  using  the  top  and 
bottom  histograms.  Figures  3  and  4  show  the  estimates  of  P^q  (m^)  arRl 
All  events  located  deeper  than  100  km  by  USCGS  were  excluded  when  obtaining  esti¬ 
mates. 

Figure  5a  shows  spectral  ratio  for  all  events  in  the  data  base,  except  those 
reported  deeper  than  100  km  by  USCGS.  Figure  5b  shows  spectral  ratio  for  only 
those  events  for  which  MSR  allows  a  decision.  These  two  figures  display  the  data 
available  to  estimate  the  probability  that  a  decision  can  be  made  using  MSR.  In  no 
case  did  an  interfering  event  destroy  the  value  of  spectral  ratio.  Figure  3  shows  the 
estimated  probability  P^^m^),  that  MSR  admits  a  decision  given  that  an  event  is  an 
earthquake,  and  no  other  event  is  interfering  with  the  short  period  data. 

Estimates  of  both  P^,q  (m^)  and  P'^q  (mb)  have  been  shown  on  Figure  3.  There 
is  no  basic  reason  why,  for  a  given  ,  these  probabilities  should  equal  each  other. 
Nevertheless  it  is  clear  that  within  the  bounds  of  experimental  error,  we  can,  in  fact, 
assume  they  are  equal.  A  nominal  curve  which  can  be  used  for  either  has  been  indica¬ 
ted  on  the  figure.  The  nominal  curve  is  the  integral  of  a  normal  distribution  with  mean 
4.  35  and  standard  deviation  0.  35.  This  curve  is  not  optimal  or  best  in  any  sense  except 
that,  subjectively,  it  appears  to  be  a  good  fit  to  the  data.  This  normal  curve  will 


henceforth  be  used  for 


and  wv- 
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Several  possible  nominal  curves  for  P  (m^)  have  been  drawn  on  Figure  4.  The 
scatter  of  experimental  data  is  quite  large.  For  this  reason  the  choice  of  nominal 
Pjg(m^)  requires  somewhat  more  comment  than  did  that  for  P^q (m^ )  or  P'^(m^). 

It  has  been  often  noted  (11)  that  on  the  average  for  any  fixed  time  interval,  the 

number  of  earthquakes  occuring  throughout  the  world  at  above  some  magnitude 

*Bmb 

tends  to  be  given  by  A10  ,  where  B  is  a  positive  parameter  on  the  order  of  unity. 

It  is  this  observation  and  the  assumption  that  earthquakes  occur  independently  accord¬ 
ing  to  a  Poisson  distribution  (12)  which  has  suggested  the  form  of  P^(m^)  used  here. 

If  A(m^)  is  the  average  number  of  earthquakes  during  a  time  interval  T  having  mag¬ 
nitudes  at  least  as  large  as  m^  then  1  -e  <V  is  the  probability  that  at  least  one 
such  event  will  occur  during  time  T.  The  several  curves  on  Figure  4  are  plots  of 
this  probability  with  A  (m^)  =  A10  ^mt).  The  large  deviation  of  experimental  points 
for  large  magnitudes  is  most  likely  a  peculiarity  of  our  very  limited  data  base. 

All  the  curves  on  Figure  4  have  been  adjusted  for  an  average  of  either  0.  04  or 
0.  08  events  having  —  5.  0  during  an  interval  T.  If  T  is  15  minutes  this  is  the 
same  as  3.  8  and  7.  6  such  events  per  day.  For  T  =  30  minutes  this  is  the  same  as 
1.  9  and  3.  8  per  day.  These  values  are  consistent  with  the  number  of  such  events 
reported  daily  by  the  USCGS.  Curves  have  been  shown  for  B  =  1. 0  and  B  =  1.  5.  These 
are  typical  B  values  one  might  find  cited  in  the  literature.  The  value  B  =  1.  5  is 
more  consistent  with  the  original  work  of  Gutenberg  (11). 
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Figure  6  is  a  summary  of  the  decision  probabilities  for  Mg  -  and  for  MSR. 

In  the  case  of  spectral  ratio  the  figure  has  been  drawn  showing  a  very  minor  effect 
of  interfering  events.  The  small  probability  of  interference  was  arbitrarily  assigned 
but  is  not  inconsistent  with  observations.  The  relatively  small  probability  of  inter¬ 
ference  for  MSR  is  due  at  least  in  part  to  the  short  time  which  must  be  free  of 
interfering  events.  For  MSR  this  time  is  usually  a  few  tens  of  seconds.  For  M  -m,  , 
the  time  may  be  15-30  minutes.  The  effect  of  interfering  events  in  the  case  of 
Mg-m^  has  been  shown  for  the  case  A  =  1. 26  x  10^  and  B  =  1.  5.  When  the  effect 
of  interfering  events  is  taken  into  consideration  we  note  that  the  probability  of  appli¬ 
cability  of  Mg-  m^  to  earthquakes  at  a  given  m^  level  is  the  same  as  that  of  MSR  for 
an  m^  value  0. 1  to  0.  2  magnitude  units  smaller. 


7 


III.  Probability  that  Mg-  m^  and  MSR  Can  Be  Applied  to  Explosions 

Let  P'  (m  )  be  the  probability  that  MSR  will  generate  a  decision  for  an  under- 
L  b 

ground  nuclear  explosion  with  body  wave  magnitude  m^.  All  of  the  limited  data  avail¬ 
able  to  estimate  P'  (m  )  is  contained  in  Figures  5A  and  5B.  P'  (m  )  has  been 
Lb  Lb 

estimated  as  follows.  The  number  of  explosions  on  Figures  5A  and  5B  in  intervals 

m,  ±0.  4  was  determined  for  m,  at  multiples  of  0. 1  magnitude  units.  Then  P'  (m,  ) 
b  b  E  b 

was  estimated  as  the  ratio  of  these  numbers.  Such  estimates  are  shown  as  the  solid 
points  on  Figure  7. 

Let  P  (m^)  be  the  probability  that  Mg  -  will  generate  a  decision  for  an 
explosion  of  magnitude  m^  .  Figure  1  contains  part  of  the  data  needed  to  estimate 
P  (m^).  That  figure  shows  all  explosions  in  our  data  base  for  which  Mg  -  m^  gen¬ 
erated  a  decision.  In  addition,  there  were  seven  explosions  for  which  long  period 
digital  data  were  available  but  which  could  not  be  identified  at  LASA.  The  LASA  body 
wave  magnitudes  of  these  events  were  5.  28,  4.  95,  4.  90,  4.  69,  4,  59,  4.  40  and  3.  86. 
The  three  unidentified  events  with  magnitudes  5.  28,  4.  59,  and  4.  40  were  the  only 

explosions  obscured  by  interfering  events.  Using  this  data  P  (m  )  has  been  estimated 

L  b 

as  the  fraction  of  all  explosions  in  the  interval  m^  ±0.  4  which  appear  on  Figure  1. 

This  estimate  has  been  evaluated  at  each  tenth  of  a  magnitude  and  is  shown  as  open 
dots  on  Figure  7. 

A  considerably  more  indirect  way  to  estimate  P^m^)  has  also  been  employed. 
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First  consider  only  explosions  without  interfering  events  and  for  which  the  surface 
waves  are  obscured  by  LP  microseismic  noise.  The  value  of  decision  probability  for 
this  set  of  events  can  be  determined  from  a  study  of  noise  with  no  events  present.  The 
M  bounds  imposed  by  noise  have  been  measured  for  a  distribution  of  hypothetical 

J 

epicenters  in  the  Sino -Soviet  Bloc  using  more  than  40  noise  samples.  Figure  8  is  a 
cumulative  histogram  of  these  bounds.  Using  this  and  the  decision  line  on  Figure  1 
we  have  calculated  the  probability  (m^)  shown  on  Figure  9.  P^  (m^)  is  just  the 
probability  that  the  M  bound  imposed  by  noise  will  locate  an  event  in  the  explosion 

O 

region.  If  the  bound  does  not  constrain  the  event  to  that  region  then  no  decision  can 
be  made.  The  probability,  P^  (m^),  that  can  be  measured  given  there  is  no  inter¬ 
fering  event  is  also  shown  on  Figure  9.  It  was  estimated  using  ratios  of  explosion 

counts  in  intervals  m,  ±  0.  4.  The  smooth  curve  for  P„  (m,  )  is  the  same  as  that  shown 

b  zb 

on  Figure  3  with  |jl  increased  from  4.  35  to  5.  45.  This  shift  of  1.  1  magnitude  units 
results  from  the  similar  difference  between  the  means  of  the  earthquake  and  explosion 
populations  shown  on  Figure  1.  If  the  probability  of  interfering  events  is  zero  then 

(m^P 

indirect  estimate  of  P^,  (m^).  It  is  very  nearly  equal  to  P^  (m^).  The  effect  of  con¬ 
sidering  interfering  events  would  be  to  shift  the  indirect  estimate  slightly  to  the  right 
in  the  same  way  that  it  shifted  the  probability  for  earthquakes  as  shown  on  Figure  6. 


inis  nas  oeen  drawn  onto  nigure  /  as  me 


PE(mb>=P2(mb)+(1  ‘P2 
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IV.  Comparison  of  Mg  -  and  MSR 

Figures  6  and  7  summarize  the  principal  quantitative  results  in  this  report.  In 
one  sense  these  figures  indicate  the  relative  value  of  the  Mg  -  and  MSR  discrimi¬ 
nants,  using  LASA  data  only,  as  a  function  of  magnitude  and  event  type.  If  an  event 
is  an  earthquake  of  magnitude  m^  it  is  clear  from  Figure  6  that  it  has  a  slightly 
higher  a  priori  probability  of  being  identified  by  MSR  than  by  Mg  -  m^  .  However,  it 
is  fair  to  say  that  the  ability  of  MSR  to  operate  at  lower  magnitudes  than  M„  -m  on 

o  D 

earthquakes  is  slight  and  amounts  to  a  shift  of  only  0. 1  -0.2  body  wave  magnitude 
units.  Roughly,  the  probabilities  in  both  cases  drop  from  1.  0  to  0.  0  as  M^  goes  from 
5.  0  to  4.  0.  The  situation  is  considerably  changed  in  the  case  of  explosions.  For 
example  it  is  clear  from  Figure  7  that  MSR  can  be  as  effective  at  —  4.  5  as 
Mg  -  can  be  at  m^  =  5.  0.  At  =4.5  MSR  will  identify  10%  of  the  explosions. 
(Presumably  no  decision  is  made  concerning  the  remaining  explosions),  whereas 
Mg  -  will  identify  virtually  no  explosions.  The  probability  for  MSR  is  similar  for 
earthquakes  and  explosions  but  for  Mg  -  it  drops  from  1. 0  to  0.  0  roughly  as 
goes  from  5.  2  to  4.  5.  The  net  impact  of  Figures  6  and  7  is  that  MSR  can  play  a 
particularly  important  role  in  discrimination  for  events  with  less  than  about  5.  0. 

Although  many  events  identified  by  MSR  are  also  identified  by  Mg  -  and  vice 
versa  the  two  criteria  can  be  quite  complementary.  This  is  true  for  large  as  well  as 
small  magnitude  events.  Consider  the  four  shallow  earthquakes  with  m^  —  4.  9 
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indicated  on  Figure  2  which  could  not  be  identified  by  Mg  -m^  due  to  interfering  events. 

A  check  of  these  events  has  shown  that  these  earthquakes  were  correctly  identified  by 

MSR.  This  is  to  be  expected  since  the  events  are  large  and  interference  of  LP  surface 

waves  will  not  generally  effect  the  short  period  P  wave  signal  to  noise  ratio.  Thus 

MSR  operates  as  successfully  as  if  there  had  been  no  LP  interference. 

A  similar  situation  holds  for  smaller  magnitude  events  as  well.  That  is,  the 

inability  to  come  to  a  decision  using  Mg  -m^  and  MSR  are  quite  independent.  This 

has  been  checked  using  earthquakes  with  —  4.  5  which  might  have  been  identified 

using  Mg  -  m^.  Figure  10  is  a  Venn  diagram,  with  areas  proportional  to  the  number 

of  events  in  the  indicated  sets,  which  shows  the  degree  of  independence  of  M  -  m, 

S  b 

and  MSR.  Both  Mg-  m^  and  MSR  identified  30%  of  the  earthquakes  considered.  If  the 
two  criteria  operate  independently  then  we  would  expect  9%  of  the  events  to  be  identi¬ 
fied  by  both  on  the  average.  In  fact  1 5%  were  identified  by  both.  Considering  the 
small  number  of  events  involved  this  is  a  strong  indication  that  both  criteria  operate 
independently  at  these  small  magnitudes.  This  independence  is  anticipated  if  it  is 
assumed  that  LP  and  SP  interference  by  noise  and  signals  is  independent  and  that 
decisions  are  withheld  whenever  noise  levels  become  too  high. 
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Fig.  2.  Surface  wave  detection  data  for  Sino -Soviet  earthquakes 
recorded  at  LASA. 
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Fig.  3.  Probability  that  Mg  -  or  MSR  can  be  applied  to 
earthquakes  with  no  interfering  events. 


16 


INTERFERENCE  PROBABILITY 
OF  UP  EARTHQUAKE  DATA 


Fig.  4.  Probability  of  interference  of  LP  earthquake  data. 
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SPECTRAL  RATIO 


Fig.  5a.  Spectral  ratio  test  results. 
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Fig.  5b.  Modified  spectral  ratio  test  results. 


19 


PROBABILITY  THAT  DISCRIMINANT  CAN 
BE  APPLIED  TO  EARTHQUAKE 


18-6-  9949 

-  BOTH  M,-mb  AND  MODIFIED  SPECTRAL  RATIO 

GIVEN  THERE  IS  NO  INTERFERING  EVENT 

- SPECTRAL  RATIO  ASSUMING  -0.01 

PROBABILITY  OF  AN  SP  INTERFERING  EVENT 

- Ms-mb  ASSUMING  POISSON  DISTRIBUTION 

OF  INTERFERING  EVENTS  WITH  AVERAGE 
NUMBER  OF  INTERFERING  EVENTS 

EQUAL  TO  1.26  x  106M,5m*> 


4r0  4,5  5.0  5.5 


LASA  mb 


Fig.  6.  Probability  that  discriminants  can  be  applied  to  an  earthquake. 
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•  DIRECT  ESTIMATE  FOR  MODIFIED  SPECTRAL  RATIO 
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7.  Probability  that  discriminants  can  be  applied  to  an  explosion. 


Fig.  8.  Apparent  surface  wave  magnitude  of  LP  noise 
(Sino -Soviet  locations  assumed). 
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Fig.  9.  Component  probabilities  for  indirect  estimate  of  probability 
that  M,,  -  m,  can  be  applied  to  an  explosion. 
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Fig.  10.  Venn  diagram  indicating  the  independence  of  M  -  m.  and 

MSR  for  earthquakes  with  m,  <4.5.  bn 
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